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ABSTRACT 

We present sub-mm photometry for 11 Hyperluminous Infrared Galaxies 
(HLIRGs, L IR > 10 13 h^L Q ) and use radiative transfer models for starbursts and 
AGN to examine the nature of the IR emission. In all the sources both a starburst 
and AGN are required to explain the total IR emission. The mean starburst fraction is 
35%, with a range spanning 80% starburst dominated to 80% AGN dominated. In all 
cases the starburst dominates at rest-frame wavelengths longwards of 50/im, with star 
formation rates > SOOMqVT . The trend of increasing AGN fraction with increasing 
IR luminosity observed in IRAS galaxies is observed to peak in the HLIRG popula- 
tion, and not increase beyond the fraction seen in the brightest ULIRGs. The AGN 
and starburst luminosities correlate, suggesting that a common physical factor, plau- 
sibly the dust masses, govern the luminosities of starbursts and AGN in HLIRGs. Our 
results suggest that the HLIRG population is comprised both of ULIRG-like galaxy 
mergers, and of young galaxies going through their maximal star formation periods 
whilst harbouring an AGN. The detection of coeval AGN and starburst activity in our 
sources implies that starburst and AGN activity, and the peak starburst and AGN 
luminosities, can be coeval in active galaxies generally. When extrapolated to high-z 
our sources have comparable observed frame sub-mm fluxes to sub-mm survey sources. 
At least some high-z sub-mm survey sources are therefore likely to be comprised of 
similar galaxy populations to those found in the HLIRG population. It is also plausible 
from these results that high-z sub-mm sources harbour heavily obscured AGN. The 
differences in X-ray and sub-mm properties between HLIRGs at z ~ 1 and sub-mm 
sources at ~ 3 implies some level of evolution between the two epochs. Either the 
mean AGN obscuration level is greater at z ~ 3 than at z ~ 1, or the fraction of 
IR-luminous sources at z ~ 3 that contain AGN is smaller than that at z ~ 1. 

Key words: infrared: galaxies - galaxies: active - galaxies: Seyfert - galaxies: star- 
burst - Quasars: general 



1 INTRODUCTION 

One of the most important results from the Infrared Astro- 
nomical Satellite (IRAS) all sky surveys was the detection of 
a new class of galaxy where the bu lk of the bolometric emis- 
sion lies in the inf rared waveband ( soifer et al. 1984 ; Banders 
& Mirabel 1996). This population, termed 'Luminous In- 



At the brightest end of the LIG population lie the Hy- 
perluminous I nfrared Galaxies (HLIR Gs), those with Lm > 
10 1 



13,0 K^Lq ( |Rowan- Robinson 2000| ). The first HLIRG to 
be fou nd (P09104+4109) was identified by Kleinmann et al. 
with a far infrared luminosity of 1.5 x 10 13 /i^ 2 Lq. 

identified F10214+4724 at 
far infrared luminosity of 



In 1991, Rowan- Robinson et al. 
z = 2.286, with an apparent 



frared Galaxies' (LIGs) , becomes the dominant extragalactic 3 x 10 ft so Lq . Later observations of this object revealed a 



popula tion at IR luminosities above 10 Lq, with a higher huge mass of molecular gas (10 h^ r [Mq ( Brown fc vanden 



space dmsity than all other classes of galaxy of comparable Bout 1991; Solomon, Downes & Radford 1992)), a Seyfert 



bolometric luminosity. 
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emissio n spectrum ([Elston ct al. 1994 ), high optical polari 
sation (Lawrence et al. 1993), and evidence for lensing with 



a mag nification of about 10 at infrared wavelengths (Gra- 
ham 



Liu 1995; Broadhurst & Lehar 1995 



al. 199£ ; Green & Rowan- Robinson 1996). These objects ap- 



Eisenhardt et 



peared to presage an entirely new class of infrared galaxy. 

The source and trigger of the IR emission in HLIRGs is 
currently the subject of considerable debate. HLIRGs may 
simply be the high luminosity tail of the ULIRG popula- 
tion, where mergers between evolved g alaxies trigger dust 



Mirabel 



shro ude d starburst and AGN activity ( Banders 
1996). There is evidence from HST observations (Farrah et 
al. 200i| that at least some HLIRGs are merging galaxies 
A second possibility is that HLIRGs m ay be very young, 
or 'primeval' galaxies. Rowan- Robinson (2000) argues that 
the majority of the emission at rest-wavelengths > 50/xm in 
HLIRGs is due to starburst activity, implying star formation 
rates > lOOOMQyr -1 . If the rest-frame far infrared and sub- 
mm emission from HLIRGs is due to star formation, then the 
star formation rates would be the highest for any objects in 
the Universe. This would strongly suggest these galaxies are 
going through their maximal star formation periods, imply- 
ing that they are galaxies in the first stages of formation. A 
final possibility is that the IR emission arises via some other 
mechanism (e.g. a transient IR luminous phase in QSO evo- 
lution not triggered by interactions), HLIRGs would then 
be an entirely different class of object. 

In this paper we study the infrared emission from 
HLIRGs using data from the optical to the sub-mm. We 
present new sub-mm data for a sample of 11 HLIRGs, and 
use radiative transfer models for starbursts and AGN, in 
conjunction with previously published IR photometry, to 
examine the power source behind the IR emission. Sample 
selection, observations and data analysis are described in 
§2. The radiative transfer models used to evaluate the IR 
emission are described in §3. Results are presented in §4, 
and notes on individual sources are given in §5. Discussion 
is presented in §6 and conclusions are summarized in §7. We 
have taken Ho = 65 km s _1 Mpc -1 , Q — 1.0 and Q.a = 0.0. 



2 THE SAMPLE 

2.1 Sample Selection and Observations 

The objects in our sample are taken fr om th e FIR selected 
sample presented by Rowan- Robinson ( 2000 ) . Unlike nearly 
all previous studies of HLIRGs, the sample in our study is 
selected in a manner independent of obscuration, inclination 
or AGN content. Together with the statistical homogeneity 
and completeness of the parent samples, our sample is there- 
fore entirely free from AGN bias and suitable for drawing 
global conclusions about the HLIRG population. 

Observations were made using the Submillimetre Com- 
mon User Bolometer Array (SCUBA, Holland et al. 1999) 
on the James Clerk Maxwell Telescope (JCMT) on Octo- 
ber 11-12 2000 and on January 8-17 2001. SCUBA contains 
two bolometer arrays, one containing 37 pixels and opti- 
mized for observations at 850/im, and the other containing 
91 pixels and optimized for observations at 450/im. In most 
circumstances, both arrays are operated simultaneously us- 
ing a dichroic. Observations were performed using SCUBA's 



photometry mode, in which data is taken using only a sin- 
gle pixel on each array. For each integration the secondary 
mirror was jiggled so that the selected bolometer in each 
array sampled a 3 x 3 grid with 2" spacing between grid po- 
sitions, centred on the source. During each integration the 
secondary mirror is chopped by 45" in azimuth with a fre- 
quency of 7Hz in order to remove sky variations. After each 
integration, the telescope is then nodded to a reference posi- 
tion 45" away in azimuth to remove hotspots in the internal 
SCUBA optics. Each object was observed for approximately 
40 minutes, depending on weather conditions. Sky opacities 
during the observations were of moderate quality, with mea- 
sured opacities at 225Ghz from the Caltech Submillimetre 
Observatory (CSO) in the range 0.075 < r 22 5 < 0.15. Cal- 
ibration observations were made of Mars or Uranus, or of 
a secondary calibrator if no primary calibrator was avail- 
able. Skydips were taken before and after each object and 
calibrator observation. 



2.2 Data Reduction 

The SCUBA User Reduction Facility (SURF) software was 
used to reduce the data for all objects. The data were first 
flatfielded and despiked. Atmospheric extinction corrections 
at 450/mi and 850/mi were derived by extrapolating from 
the CSO T225 extinction valu es fol lowing the prescription of 
Archibald, Wagg & Jenness (2000). These extrapolated val- 



ues were checked for consistency against the observed 450/im 
and 850/mi extinctions from the skydips. Residual sky gra- 
dients not removed by nodding and chopping were removed 
by averaging over all the bolometers in each array, and sub- 
tracting this value from the measured source flux. Individ- 
ual integrations for each source were then concatenated into 
a single exposure. Each concatenated dataset was checked 
for internal consistency using a Kolmogorov-Smirnov (K-S) 
test. Finally, flux calibration for each source was carried out 
using the FLUXES package together with the calibrators 
listed in Table |l|. In those cases where the calibrator has 
a larger spatial extent than the 9 point photometry mode 
jigglemap the flux of the calibrator in a single photometry 
jigglemap was used to calibrate the sources, rather than the 
total flux of the calibrator. 



3 INFRARED EMISSION MODELS 
3.1 Starburst Models 

To model the IR emission due to starburst activity we 
used the starburs t mo dels of Efstathiou, Rowan- Robinson 
& Siebenmorgen (2000). These models consider an ensemble 
of evolving HII regions containing hot young stars, embed- 
ded within Giant Molecular Clouds (GMCs) of gas and dust. 
The composition of the dus t is giv en by the dust grain model 
of Siebenmorgen & Krugel (1992), which includes Polycyclic 



Aromatic Hydrocarbons. The stellar populations within the 
GMCs evolve accor ding t o the stellar synthesis codes of 
Bruzual & Chariot ( 1993 ). The star formation rate is as- 
sumed to decay exponentially with an e-folding timescale of 
2 x 10 7 years. The models vary in starburst age from zero 
years up to 7.2 x 10 7 years, with 11 discrete values. 
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Table 1. Hyperluminous Galaxies observed by SCUBA 



Name RA (J2000) Dec z Type 1 Obs. Date Calibrator 

hh:mm:ss ° ' " 



IRAS F00235+1024 


00 


26 


06.7 


10 


11 


27.6 


0.58 


nl 


October 


2000 


Uranus 


IRAS 07380-2342 


07 


10 


09.8 


-23 


49 


57.9 


0.29 


ul 


January 


2001 


oh 231 


IRAS F 10026+4949 


10 


05 


52.5 


49 


34 


47.8 


1.12 


Syl 


January 


2001 


oh 231 


IRAS F12509+3122 


12 


53 


17.6 


31 


05 


50.5 


0.78 


QSO 


January 


2001 


Mars 


IRAS 13279+3401 


13 


30 


15.3 


33 


46 


28.7 


0.36 


QSO 


January 


2001 


Mars 


IRAS 14026+4341 


11 


04 


38.8 


43 


27 


07.2 


0.32 


Syl 


January 


2001 


Mars 


IRAS F14218+3845 


11 


23 


55.5 


38 


31 


51.3 


1.21 


QSO 


January 


2001 


Mars 


IRAS F16124+3241 


16 


14 


22.1 


32 


34 


03.7 


0.71 


nl 


January 


2001 


Mars 


ELAISP90 J164010+410502 


16 


10 


10.2 


41 


05 


22.1 


1.10 


QSO 


January 


2001 


Mars 


IRAS 18216+6418 


18 


21 


57.3 


61 


20 


36.4 


0.30 


Syl 


October 


2000 


Uranus 


IRAS F23569-0341 


23 


59 


33.6 


-03 


25 


12.8 


0.59 


nl 


October 


2000 


Uranus 



Positions and redshifts are taken from the NASA Extragalactic Database. 1 Spectral Type, taken from Rowan- 
Robinson 2000 & NED. 'nl' - Narrow Line object, 'Syl' - Seyfert 1. 



3.2 AGN Torus Models 

To model the IR emission due to AGN, we used t he AGN 
models of Efstathiou &i Rowan- Robinson (1995). These 
models incorporate accurate solutions to the axially sym- 
metric radiative-transfer problem in dust clouds to model 
the IR emission from dust in active galactic nuclei. Dust 
compositio n is g iven by the multigrain dust model of Rowan- 
Robinson (1992). We have used the thick tapered disk mod- 
els following an r _1 density distribution, as this subset of 
models has been found to be most successful in satisfying the 
observational constraints of AGN. The AGN models vary in 
torus opening angle from 0° to 90°, with 15 discrete values. 



4 RESULTS 

4.1 Spectral Energy Distributions 

We combined our measured sub-mm fluxes and 3<r upper 
limits with optical and IR data from the literature to fit 
Spectral Energy Distributions (SEDs) for each object. For 
each source we obtained IRAS flux es and 3a upper li mits 
from the Faint Source Catalogues (Moshir et al. 1992), or 
by using the SCANPI v5.0 software. Where available, we ob- 
tained B and R band magnitudes from the APM catalogues, 
and near-IR fluxes from the 2MASS catalogues. Additional 
photometry for each source is described in section 6. The 
compiled fluxes for each object are given in Table 

Emission from unobscured population II and III stars 
is not included in either the starburst or AGN models, and 
care must be taken in accou nting for this. For example, a 
recent HST study of ULIRGs ( [Farrah et al. 200 lj ) found that 
the optical emission was in most cases dominated by old stel- 
lar populations and not by light from a starburst or AGN. 
To avoid the uncertainties in assuming an arbitrary SED 
for an evolved stellar population, we have assumed that all 
emission shortwards of 4^im in the rest frame of the objects 
contains a significant but unquantified contribution from old 
stellar populations. As this contribution could lie between 
0% and 100%, any measured flux at a rest-frame wavelength 
shorter than 4fim is treated as a 3er upper limit in the fit- 
ting, unless the contribution from the host galaxy has been 
previously removed. 



Goodness of fit was examined by using the reduced \ 2 
statistic. Fits for all sources were good, with Xbest — 1-5 m 
all cases. To determine the errors on the luminosities and 
model parameters we explored x 2 space between xlest an d 
Xbest + 2. The SEDs for those sources with sufficient IR 
photometry to constrain the shape of the SED are given 
in Figure |jj The total IR luminosities for each source, the 
starburst and AGN components, and the best fit model pa- 
rameters are given in Table ^. 

There exist models for the IR emission from dusty QSOs 
( Rowan- Robinson 1995) where the dust can extend several 



hundred parsecs from the nucleus inst ead of the several tens 



of parsecs found in other AGN models ( Efstathiou fc Rowan- 



Robinson 1995). The dust at these large radii will have a 
temperature of 20 - 70 Kelvin, as opposed to the dust several 
tens of parsecs from the AGN, which will have a tempera- 
ture of a few hundred Kelvin. If extended dust around an 
AGN were present in any of our sample the sub-mm emission 
from this dust could be mistakenly ascribed to star forma- 
tion, thus overestimating the star formation r ates. To test 
this, we fitted the models of Rowan- Robinson ( |l99.j ) to the 
objects in our sample, both on their own and in combina- 
tion with the starburst models described previously. In all 
cases the extended AGN dust models were clearly rejected 
by the fitting. We therefore conclude that sub-mm emission 
from extended dust surrounding an AGN is not viable for 
the objects in our sample. 



4.2 Star Formation Rates 

Estimating obscured star formation rates from IR data is 
based on silicate and graphite dust grains absorbing the op- 
tical and UV light from young stars and reradiating in the 
IR and sub-mm. A recent estimate for deriving star forma- 
tion rates from IR lu minosities has been made by Rowan- 
Robinson et al. (|1997|): 



M*,aii= 2.6 x 10 -— - 



(1) 



where M*, a ii is the rate of formation of stars, e is the frac- 
tion of optical/UV light from the starburst that is absorbed 
by the dust and re-emitted in the IR, and Leo is the 60/im 



© 2001 RAS, MNRAS 000, ??-?? 



4 D. Farrah et al. 



Table 2. Compiled Photometry from the optical to the sub-mm for the HLIRGs in our sample 



Name 


Si 


Sj 


Sh 


Sk 


^6.7 


Sl2 


S25 




S100 


5450 


■5850 


F00235+1024 


0.08 








0.92 ± 0.38 


< 173 


< 193 


428 ± 55.6 


< 938 


< 967 


< 6.38 


07380-2342 6 


1.11 










484 ± 28 


801 + 80 


1 1 7(1 + 88 


3546 + 280 


< 1275 


26.8 ± 4.2 


F10026+4949 


0.04 










< 86 


177 ± 35 


266 ± 30 


310 ±80 


< 386 


< 7.91 


F12509+3122 c 


0.86 


0.96 


0.81 


1.13 




< 106 


103 ± 26 


218 ± 44 


309 ± 82 


< 333 


< 9.23 


13279±3401 d 












< 94 


< 126 


1182 ± 77 


1196 ± 165 


< 274 


< 9.22 


14026+4341 e 




3.36 


4.76 


8.19 




118 ±27.1 


285 ± 29 


622 ± 57.0 


994 ± 239 


< 94 


< 7.53 


F14218+3845 / 


0.05 








0.79 ± 0.26 


< 96.9 


< 74.9 


< 565 


< 2100 


< 253 


< 8.55 


F16124+32419 








0.15 




< 65 


< 55 


174 ± 35 


290 ± 90 


< 103 


8.47 ± 2.2 


E J1640+41' 1 




0.44 


0.63 


0.56 


4.99 ± 0.66 


< 102 


< 57 


< 81 


< 336 


< 82 


< 6.97 


18216+6418* 




8.8 


13.0 


22.1 


109 ± 11 


211 ± 14 


395 ± 45 


1128 ± 113 


2159 ± 252 


< 308 


14.8 ±2.6 



All fluxes are given in mjy. Unless stated in the text, all fluxes at rest-frame wavelengths shortward of Afim are treated as upper limits 
in the SED fitting, due to the unquantified contribution from the host galaxy. Sources for additional photometry are described in §6. 
Additional photometry: (a) Si 5 = 6.752 ± 2.143, S 90 = 477.9 ± 147.9, Si 80 = 804.7 ± 393.3 (b) S B = 0.062, S R = 0.346 (c) S B = 1.41, 

S v = 0.751, S R = 0.72 (d) S B = 0.50, S R = 2.38 (e) S B = 3.23, S R = 5.48, S4.0 = 16, S7.0 = 36 ± 4 (f) S B = 0.06, S R = 0.114, 
S15 = 3.228 ± 1.036, S 90 = 163.5 ± 61, Si 80 =< 1765 (g) S B = 0.014, S R = 0.079 (h) S15 = 7.75 ± 0.72, S 90 = 72 ± 23 (i) Su = 9.04, 
S s = 8.12, S v = 7.55, S L = 60.5, S9.63 = 143 ± 14.3, Snoo =< 47, S1250 =< 3.7 



starburst luminosity. The factor (f> incorporates the correc- 
tion between a Salpeter IMF and the true IMF (xl.O for 
a Salpeter IMF, x3.3 for a Miller-Scalo IMF), and also a 
correction for the assumed upper and lower stellar mass 
bounds for the stars forming in the starburst (xl.O if the 
mass range is 0.1 < Mq < 100, x 0.323 if the mass range is 
assumed to form only OBA type stars, i.e. 1.6 < Mq < 100). 
We have assumed that all of the optical/UV light from the 
starburst is absorbed (e ~ 1.0), and that the IMF of the 
starburst is a Salpeter IMF forming stars across the mass 
range 0.1 < Mq < 100, (<j> = 1.0). The calculated star for- 
mation rates for our sample using Equation ^ are given in 
Table |. 



4.3 Dust Masses and Temperatures 

Dust temperatures cannot be directly determined from the 
starburst or AGN models described in sections 3.1 and 3.2, 
as these models do not assume an isothermal temperature 
distribution. Monolithic dust temperatures, although un- 
physical simplifications in AGN, can serve as a useful com- 
parison with previous work. Dust temperature estimates 
were calculated by fitting an optically thin greybody func- 
tion of the form: 



S u = v^B{v,T d u 



to the FIR SED over the wavelength range 200 - 1000/im. 
In Equation ^| S v is the source flux at a frequency v, Tdust 
is the dust temperature, B(v, Tdust) is the Planck function, 
and 13 is the frequency dependence of the grain emissivity. 

Gas and dust masses can be estimated directly from the 
starburst models. The total mass of Giant Molecular Clouds 
(CMC's, which include both gas and dust) is straightforward 
to estimate, as stars are assumed to form from the CMC's. 
The mass of stars at time t* can be written: 



M a (f) = rjMoMc(f) = / Ms (t)dt 



(3) 



where MGMc(t*) is the mass of Giant Molecular Clouds at 
time t* , Ms {t) is the star formation rate and rj is the effi- 
ciency of conversion of CMC's to stars. Approximating the 



star formation rate as an exponential decay with a charac- 
teristic e-folding time of r = 20Myr (ERS 00), then it is 
trivial to show that: 



M G Mc(t*) 



r Ms (t*) 



(4) 



where Mgmc is given in multiples of 1O 6 M as r is given m 
multiples of 10 6 years. This however gives the total mass of 
GMCs in the system at the current age of the starburst, 
without accounting for the effects of supernovae winds, 
which will remove some fraction of gas and dust from the 
system either via photodestruction of dust grains or by blow- 
ing the gas and dust out of the starburst regions. Estimates 
using Equation ^ may therefore significantly overestimate 
the gas and dust masses in the system, unless the starburst 
is very young. For the purposes of this analysis it is bet- 
ter to estimate dust masses directly from the sub-mm SED. 
An approach for estimating dust and gas masses in galax- 
ies using sub-mm data is to assume the system is optically 
thin at these wave lengths, and to follow the prescription of 
Hildebrand (^98^): 



M d 



1 + Z K,(v r )B(v r , Tdus 



(5) 



where v a and v r are the observed and rest frame frequen- 
cies respectively, S Vo is the flux in the observed frame, 
B(v r , Tdust) is the Planck function in the rest frame and 
Tdust is the dust temperature. The gas mass is then ob- 
tained by assuming a fixed gas to dust ratio. For the most 
extreme IRAS galaxies, th e best current estimate of the ga s 
to dust ratio is 540 ± 290 flSanders, Scoville fc Soifer 199l[). 



For comparison, the gas to dust ratio in spira l galaxies is 
thought to b e ~ 500 (Devereux k Young 1990| ), and ~ 700 
in ellipticals (Wiklind & Henkel 1995). The mass absorption 
coefficient in the rest frame, n r , is taken to be: 



(6) 



in units of m 2 kg 1 and the luminosity distance is assumed 
to be of the form: 



D L = 



H q 2 



(q Q z + (q Q ~ l)[(2q z + I)* - 1]) 



(7) 
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Dust masses and temperatures for our sample calculated 
using Equations ^ and |B| are given in Table [| 



4.4 Individual Sources 

44.1 IRAS F00235+1024 

This narrow line object is not detected in the X-ray to a limit 
of Lx/Lboi = 2.3 x 10 -4 , indicating either atypically weak 
X-ray emission or an obscuring colum n of Nh > 10 23 cm~ 2 
( [Wilman et al. 1998j ). HST imaging jFarrah et al. 2002j ) 
shows an i nteracting object, w ith no optical QSO. ISO ob- 
servations (Verma et al. 2002) imply that the IR emission 
is dominated (> 60%) by a starburst, but that an AGN 
contribution is required to explain the MIR emission. 

The best fit SED model for this object is given in Fig- 
urejll Additional / band data was taken from Farrah et al. 
(p002|) , and additional ISO fluxes were compiled from Verma 
et al. (2002). The SED of this object is best explained by 



a combination of starburst and AGN components. The ISO 
detections at 6.7pm and 15pm combined with the sub-mm 
upper limits make it impossible to fit this source with a 
pure starburst. We find that the AGN and starburst com- 
ponents contribute comparably to the total IR emission. The 
starburst component dominates at wavelengths > 60pm. 
The best fit model parameters are a dust temperature of 
35 Kelvin, a dust mass of 1.95 x 10 8 A/q, a star formation 
rate of ~ lOOOA/Qyr^ 1 , and a line of sight angle to the AGN 
torus of < 17.2°. The strong mid-IR absorption feature at 
15pm predicted by the combined SED fit results from view- 
ing an AGN dust torus almost edge on through intervening 
absorption by cooler dust. 



44.2 IRAS 07380-2342 

The best fit SED model for this object is given in Figure |l| 
Additional / band data was taken from archival HST obser- 
vations. This narrow line HLIRG is best fitted by a combi- 
nation of starburst and AGN components where the AGN 
component dominates overall (~ 65%), but where the star- 
burst dominates at wavelengths longward of 45pm. The IR 
data is best fit by a young (6.6 - 16Myr) starburst with a 
star formation rate of 1250A/Qyr _1 , and an AGN component 
with a line of sight angle of < 20°. The dust temperature 
and mass are 31K and 6.6 x 10 s Mq respectively. 



4.4.3 IRAS F10026+4949 



HST imaging of this source (Farrah et al. 2002) shows mod- 
erate morphological disturbance and multiple very close 
companions. The object is thus assumed to be in ongoing 
interactions. We compil ed ad ditional / band data for this ob- 
ject from Farrah et al. (2002). With only three fluxes, there 



exists insufficient IR photometry to model the IR SED. The 
fluxes combined with the upper limits do however allow con- 
straints to be drawn on both the total IR luminosity, and 
the starburst and AGN components. The best fit predicts a 
total IR luminosity of ~ 6.5 x 10 13 L Q , with the AGN pro- 
viding ~ 80% of the total IR emission. The star formation 
rate is extremely high, at ~ 2OOOA/0yr _1 . 



444 IRAS F12509+3122 

This QSO shows no signs of i nteraction, and has an undis- 
turbed bright elliptical host (Farrah et al. 2002). We ob 
tained additional / band data from Farrah et al. (2002) 



which is treated as a flux rather than an upper limit in 
the SED fitting as the host galaxy magnitude is known. 
The best fit to the IR SED predicts a total IR luminos- 
ity of ~ 1.82 x 10 13 L©, where the AGN component pro- 
vides the majority (~ 60%) of the emission over the whole 
SED but where the starburst dominates at wavelengths 
> 50pm. The dust temperature and emissivity are 39K and 
1.94 respectively, with a predicted star formation rate of 
~ lOSOMoyr" 1 . An SED fit is plotted in Figure §. The 
model slightly underpredicts the I band flux. This is likely 
due to treating the other optical fluxes as upper limits. As 
this object is a QSO the optical emission will arise mostly 
from the AGN, and ind eed the host galax y contribution in 
the / band is very small (Farrah et al. 2002). For consistency 
with the other sources we have treated the other optical 
fluxes as upper limits, as trial fits show that fitting to the 
optical points as fluxes makes no difference to the derived 
IR luminosities. 



44.5 IRAS 13279+3401 

This broad line object has only two measured fluxes, so con- 
straints can only be drawn on the total IR luminosity, and 
the best fit AGN and starburst components. The best fit pre- 
dicts a total IR luminosity of ~ 8.7 x 1O 12 L with the AGN 
component providing the majority (~ 70%) of the emission. 



4.4.6 IRAS 14026+4341 

This r adio quiet BAL QSO was discovered by Low et al. 

nil- 

HST imaging shows that the host is a luminous el- 



liptica l, and that the host is in volved in ongoing interac- 
tions (Hjrtchjn^s_^^foms_1^995). Optical spectropolarime- 
try (Hines et al. 2001) shows that the polarized flux density 
is dependent on wavelength, indicating that the polarization 
is produced by dust scattering. 

The best fit SED for this object is presented in Figure |l| 
For this object we obtained additional unpublished ISO pho- 
tometry (Belinda Wilkes, priv. comm.). The best fit model 
predicts a total IR luminosity of 7.9 x 1O 12 L with the AGN 
component providing ~ 60% of the IR emission. The derived 
dust temperature and mass are 35K and 1.7 x 10 8 Mq respec- 
tively, with a star formation rate of ~ 6OOA/0yr _1 . 



4.4.7 IRAS F14218+3845 

This QSO possesses a small very luminous host with no signs 
of ongoing interact ions and an undete rmined morphology 
(jFarrah et al. 2002|). ISO observations ([Verma et al. 20o|) 



suggest that the IR emission is best explained as being pre- 
dominantly (74%) starburst in origin, but with the AGN 
component dominating at wavelengths < 30pm. 

A dditio nal ISO photometry was compiled from Verma 
et al. ( 2002 ), and additional / band data from Farrah et al. 
( [2002| ). The / band data is treated as a flux rather than 
an upper limit as the host galaxy magnitude is known. The 
best fit predicts a total IR luminosity of 1.15 x 1O 13 L0. The 
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Figure 1. Best fit combined models in the observed frame, and 
F12509+3122, 14026+4341, F14218+3845 & 18216+6418 
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AGN and Starburst components, for F00235+1024, 07380-2342, 
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Figure 1 - continued 



starburst component dominates, providing 80% of the total 
IR emission. All the model fits are consistent with a star- 
burst age of < 57Myr, and a dust mass and temperature 
of 1.07 x 10 s Mq and 45K respectively. The predicted star 
formation rate is extremely high, at 2000Moyr -1 . The best 
fit SED is presented in Figure [l]. The starburst component 
dominates at rest-frame wavelengths > 45/xm, with the AGN 
component dominating in the optical. 



44.8 IRAS F 16124+3241 

This object shows no detectable polarized broad lines in the 
optical, suggesting that the IR emission may be power ed by 



a young st arburst rather than a dust shrouded AGN (Tran 



et al. 199E). With only three fluxes, it is not possible to draw 
constraints on all model parameters. The best fit total IR 
luminosity is ~ 1.1 x 10 13 Lq, with the starburst providing 
~ 60% of the IR emission. An AGN component is however 
required to explain the total IR emission over 1 — lOOO^im. 

44.9 ELAIS J1640+41 



This radio quiet QSO was dis covered (Morel et al. 2001 ) 
as part of the ELAIS survey ( piiver et al. 2000p. Unlike 



many other HLIRGs this object is detected in the soft X- 
ray, suggesting that this object m ay represent the face on 
analogues of narrow-line HLIRGs (Morel et al. 2001). 

A SED fit f or thi s object has been previously presented 
by Morel et al. (2001). Here we refit the SED with improved 



estimates for the IRAS upper limits using the SCANPI soft- 
ware. We derive a total IR luminosity of 7.9 x 10 12 Lq with 
the AGN providing ~ 70% of the IR emission. The star for- 
mation rate is 55OM0yr _1 . These results are in agreement 
with Morel et al. 



44.10 IRAS 18216+6418 

This broad line object posesses some enigmatic properties. 
Radio observations classify it as radio quiet, although the 
richness of the environm ent (Abell class ~ 2) is more typical 
of Radio Loud objects (Lacy, Rawlings & Hill 1992). Many 
other properties of the system are more typically found in 
radio loud objects The host galaxy, with a magnitude of 
Mv = —23.2, is a large, luminous elliptical with no signs 
of ongoing interaction. The QSO nucleus is slightly offset 
from the centre of the host galaxy. Radio mapping o f this 
sources show a core and 



pair of antip odean jets (Lacy. 
Rawlings fc Hill 1992) ; [Blundell et al. 1996; ), a radio structure 



more typically found in RLQs. The radio structure implies 
that a central supermassive black hole is respons ible for the 
compact radio emission, rather than a starburst ( Blundell et 



al. 1996). Off- nuclear optical spectroscopy shows a spectrum 



more typically found in radio loud objects (Fried 1998) 



In addition to our observations of the QSO we also ob- 
served with SCUBA the other object within the IRAS er- 
ror ellipse and found no detection, further supporting the 
identification of the QSO as the source of the IRAS emis- 
sion. Additional optical photometry for this object was taken 
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Table 3. Hypcrluminous Galaxies: Luminosities and Model Parameters 



Name 



j Tot 
L IR 



rSb 
L IR 

Lq 



t AGN 
L IR 



Age a 
Myr 



SFR 6 
Meyr" 1 



Mj 

a 

Mq 



kelvin 



0* 



F00235+1024 

07380-2342 

F10026+4949 

F12509+3122 

13279+3401 

14026+4341 

F14218+3845 

F16124+3241 

E J1640+41 

18216+6418 



13 04+ 08 
13 34+0.01 

11 O1+0-03 
id - 81 -n.l3 

13.261^? 
12.881°;^ 

12 90+ - 03 

13 06+ 02 

13 02+ 02 
ro.uz_ 14 

12 90+ 01 
1+0.01 
0.01 



13.14" 



12 74+0-06 
12.88^ 

12 og+0.30 

12 47+ - 14 
iz - 4 <_0.05 

12.96lg-« 
lz -' a -o.io 

12 oq+0.09 
.+0.02 



12.74 n 



0.02 



12 73+ - 11 

13 16+ 001 

,o yq+0.03 

io./o_ 13 

13 02+ 02 

19 71 +0.05 

iz - (1 -0.18 

12 70+ - 03 
iZ -' U -0.13 

12 35+0- 10 
i4.oo_ Q2s 

12 63+ - 10 

ir, yo + 0.01 

,+0.01 
0.01 



12.91" 1 



6.6 - 16 



< 57 

< 57 



26 - 37 



1000 ± 300 
1250 ± 300 
2000 ± 600 
1050 ± 250 
500 ± 300 
600 ± 150 
2000 ± 500 
1100 ± 300 
550 ± 200 
1100 ± 200 



< 17.2 
< 20 



> 26 



> 46 



8.29 
8.79_ 
< 8.7 



+0.11 
0.55 
+0.10 
0.21 



8.20 



+0.08 



8.24 
8.03 
8.54 



f 0.05 
-0.50 
f0.50 
-0.13 
+0.1 
0.4 

< 8.2 



,.51 



+0.03 
-0.02 



35 + 5 1.92 + 0.1 

31 + 3 1.97 + 0.05 

39 + 7 1.94 + 0.07 

35 + 5 1.94 + 0.02 

45 + 10 1.97 + 0.02 

32 + 5 1.94 + 0.05 

31 + 4 1.94 + 0.01 



Luminosities are the logarithm of the 1-1000/xm luminosities obtained from the best fit combined starburst/AGN model, in units 
of bolometric solar luminosities. "Starburst age 6 Star formation rate c Viewing angle of the AGN dust torus d Dust mass e Dust 
temperature ^Emissivity index 



from K olma n et al. (1993), additional ISO fluxes from Clavel 
et al. (2000), and additional sub-mm data from Andreani, 
Franceschini & Granato (199£). The comprehensive photom- 
etry for this object allows excellent constraints on the best 
fit model parameters. The best fit total IR luminosity is 
1.4 x 10 13 L Q , with the AGN component providing ~ 60% 
of the IR emission. The best fit starburst model parameters 
are a starburst age of between 26Myr and 37Myr, a dust 
mass of 3.2 x 1O 8 M0, and a dust temperature of 31K. The 
inferred star formation rate is llOOMoyr -1 . 



4.4.11 IRAS F 23 5 69- 0341 



Previous ISO observations of this object ( Verma et al. 2002) 
did not confirm the IRAS 60/im detection, and did not de- 
tect the object in any other band. As we did not detect this 
source in the sub-mm, its reality remains to be determined. 
With only one IR flux, it is impossible to draw meaning- 
ful constraints on either luminosities or model parameters, 
we therefore only present a compilation of IR data for this 
object. Worth noting however is that all the model fits are 
consistent with a total IR luminosity of lO 13 ' 5 ^© or less. 



5 DISCUSSION 

5.1 Starburst & AGN Properties 

The contributions from starburst and AGN activity to the 
total I R emission in the most luminous IRAS galaxies has 



been discussed by several authors. Mid-IR spectroscopy 
( penzel et al. 1998] ; [Rigopoulou et al. 199S| ) found that the 
IR emi ssion in most ULIRGs (~ 80%) was powered by star- 



Table 4. Additional HLIRGs from the literature with calculated 
starburst and AGN fractions 



Name 




z 


TSb 
IR 

Lq 


t AGN 
h IR 

Lq 


rTot 
h IR 

Lq 


SMMJ02399-0136 


2.803 


12.89 


12.47 


13.03 


IRAS F10214- 


-4724 


2.286 


13.10 


13.27 


13.49 


PG1148+549 




0.969 


13.27 


13.51 


13.71 


IRAS F12514- 


rl027 


0.30 


12.67 


12.63 


12.95 


H1413+117 




2.546 


12.71 


13.38 


13.46 


IRAS F14481- 


-4454 


0.66 


12.94 


13.03 


13.29 


IRAS F15307- 


h3252 


0.93 


12.96 


13.28 


13.45 



Luminosities are given for Ho = 65 km s Mpc , Qq = 1.0. 

The luminosities of SMMJ02399-0136, F10214+4724 and 
H1413+117 have been corrected for the effects of gravitational 
lensing. Taken from Rowan-Robinson (2000) & Verma et al. 
(2002). 



below 10 12 ' 4 L Q , most ULIRGs were starburst dominated, 
with the starburst component contributing around 85% to 
the IR emission. At IR luminosities above 1O 12 ' 4 L0 the AGN 
contribution was much higher, contributing at least 50% of 
the IR emission. 

Studies of samples of HLIRGs to look for contribu- 
tions from starbursts and AGN to the total IR emission 
have yielded interesting results. Studies of small samples 
of HLIRGs produced conflicting results, with some authors 
suggesting that the IR em ission in HLIRGs arises predom- 



inant ly from a starburst (Frayer et al. 1995; Frayer et al 



199£ ) with star formation rates of the order ltrMoyr 
and other aut hors suggesting that HLIRGs are powered by 



a dusty AGN ( Granato, Dancse & Franceschini 1996 



Evans 



bursts, although at least half of their samples showed evi- et al. 19*98] ; Yuri fc Scoville 1998| ). There is also evidence from 



dence for both starburst and AGN activity. There was no 
detected trend for the AGN-like systems to reside in the 
more compact (and hence more advanced merger) systems. 
The fraction of sources powered by an AGN has been shown 
by several authors to increase with increasing total IR lu 



s pectropolarimetry for a buried AGN in several HLIRGs 
( |Hines et al. 1995| ; jGoodrich et al. 1996j) Previous SEP mod - 
elling of HLIRGs (|Rowan-Robinson 2000[ |Vcrma et al. 2002;) 



minosity (Veilleux et al. 1995; Shier, Rieke & Rieke 1996 



Veilleux, Sanders fc Kim 1999 ) . Recent ISO observations of 



ULIRGs (Lran et al. 2001) found that, at IR luminosities 



found that most required an AGN and a starburst to power 
the total IR emission. Overall, about half were AGN domi- 
nated. These results were however based on samples either 
biased toward AGN or lacking the tight constraints on star- 
burst luminosities given by sub-mm data. 
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The total IR luminosities of the 10 reliably detected ob- 
jects in our sample, together with the starburst and AGN 
fractions, are presented in Table ^. In all ten objects both 
a starburst and an AGN are required to explain the total 
IR emission; in all cases a 'pure' starburst or AGN model is 
unable to explain the emission over 1 — 1000/im. The mean 
starburst fraction in our sample is ~ 35%. In all ten objects 
both the AGN and starburst components supply at least 
20% of the total IR luminosity, any of the objects would 
be classified as at least a ULIRG with only the starburst 
or AGN component present. Both starburst and AGN ac- 
tivity in HLIRGs are therefore central in interpreting their 
properties. 

Figure |s| shows a plot of bolometric AGN luminosity 
plotted against bolometric starburst luminosity for our sam- 
ple. Also included in this figure are those HLIRGs from 
Rowan-Robinson (2000) and Verma et al. (2002) that have 
constrained starburst and AGN fractions. These additional 
HLIRGs are listed in Table It should be noted that only 
two of these objects (SMMJ02399-0136 & F10214+4724) are 
selected in a manner unbiased toward AGN. Three of the 
objects in Table | are gravitationally lensed (SMMJ02399- 
0136, F10214+4724 & H1413+117), and the effects of lens- 
ing must be corrected for to allow comparisons with unlensed 
objects. For SMMJ02399-0136 we have used the lensed 
s tarbu rst and AGN luminosities given by Row an- Robinson 



(200C) and applied a lensing correction of 2.5 (Ivison et al 
1998;). for F10214+4724 we have taken the unlensed star- 
burst and AGN luminosities from Green & Rowan- Robinson 



( [I996D . For H1413+117 (the 'Cloverleaf') we have taken the 
lensed sta rburs t and AGN luminosities given by Rowan- 
Robin son ( 200C ) and a pplied a lensing correction of a factor 
of 10 ( |Yun et al. 199^ . 

It can be seen from Figure ^ that there is a wide spread 
in fractional power sources across the sample, ranging from 
80% starburst dominated to 80% AGN dominated. This 
i s con sistent with a similar plot given in Rowan- Robinson 
( 2000 ) and confirms that a wide spread in fractional star- 
burst luminosities exists over at least 5 orders of magnitude 
in total IR luminosity in the IR galaxy population. The mea- 
sured star formation rates are all greater than SOOMQyr" 1 , 
and in seven cases exceed lOOOM0yr _1 . For the six objects 
in the sample where the SED shape is constrained, the IR 
emissio n longward of 50/im in the rest-frame is in all cases 



domina ted by star formation. 



Further deductions can be made from Figure |2j. With 
only 5 systems in which the starburst contributes more to 
the total IR luminosity than the AGN it is impossible to 
draw conclusions on any correlation between starburst and 
AGN luminosity. Amongst the 12 AGN dominated systems 
however there is a strong correlation between the starburst 
and AGN luminosities. Furth er comparison with figure 19 
from Rowan- Robinson (2000) shows that this correlation 
extends over 5 orders of magnitude in total IR luminosity. 
Considering only the AGN dominated HLIRGs in Figure bl 
the Spearman rank correlation coefficient between the AGN 
luminosity and the starburst luminosity is 0.95. When the 
additional AGN dominated systems from Table ^ are in- 
cluded the correlation coefficient becomes 0.98. 

This is an interesting result for understanding AGN 
and starburst activity in IRAS galaxies at all luminosities, 
as it implies that the luminosities of the AGN and star- 



burst components are influenced by common physical fac- 
tors. On prima facie grounds this may seem counterintuitive. 
Recent studies of starburst an d QSO triggering in merg- 
ers between gas rich galaxies ([Barnes fc Hernquist 1996 



Taniguchi, Ikeuchi & Shioya 1999; Bekki 2001) suggest a 



complex pattern of starburst activity. Depending on the im- 
pact parameters of the merger, the number of merger pro- 
genitors and their morphologies, there may be multiple star- 
burst events triggered throughout the lifetime of the merger. 
AGN activity conversely, is likely to follow a simpler pattern. 
If a > 10 s Mq black hole is not located in either of the merger 
progenitors then any central AGN will not reach QSO level 
luminosities until at least 10 s years after the merger has 
commenced ( Taniguchi, Ikeuchi fc Shioya 199S ), leading to 
a potential lag between any single starburst and AGN event. 
It may thus be expected that there should be at best a weak 
correlation between starburst or AGN luminosities against 
total IR luminosities in interaction driven IRAS galaxies. 
Our results however suggest that there is an underlying fac- 
tor linking the luminosities of the starburst and AGN com- 
ponents. As the physical mechanisms behind star formation 
and AGN activity are different, the most likely candidate is 
the available quantity of fuel in the circumnuclear regions 
of the system. The same reasoning applies if we consider 
HLIRGs as young galaxies in the process of formation. 

The fraction of the IR luminosity due to AGN activity 
plotted against the total IR luminosity is given in Figure 
|§| Also plotted are the additional HLIRGs from Table It 
can be seen that there is no clear trend of increasing AGN 
fraction with increasing IR luminosity amongst the HLIRG 
population, irrespective of whether biased or unbiased sam- 
ples are selected. We therefore find that the trend observed 
in the ULIRG population of increasing AGN fraction with 
increasing IR luminosity peaks in the HLIRG population, 
and that the AGN fraction does not continue to increase at 
IR luminosities greater than W 13 Lq. 

This study has demonstrated that, to accurately gauge 
both the total IR luminosity in the most luminous IRAS 
sources, and the contributions from starburst and AGN com- 
ponents, it is necessary to have fluxes and upper limits from 
the near-IR to the sub-mm. Two objects in our sample 
(F00235+1024, F14218+3845) have different IR luminosi- 
ties, and starburst and AGN fractions, than those derived 



by p revious authors ( Rowan- Robinson 200C; Verma et al 



2002) without sub-mm fluxes or upper limits. 



5.2 Evolutionary Models for HLIRGs 

A prevalent picture for the evolution of th e most luminous 
IRAS galaxies is the Sanders et al. (1988) picture. Under 



this picture (hereafter the S88 picture), ULIRGs are the 
dust shrouded precursors to optical QSOs. Interactions and 
mergers between gas rich spirals transport gas to the central 
regions of the galaxies. This large central gas concentration 
triggers starburst activity, and in the latter stages of the 
merger commences the fuelling of a central supermassive 
black hole. In the last stages the dust screen shrouding the 
black hole is blown away and the ULIRG evolves into an op- 
t ical QSO. A more recent evolutionary scenario for ULIRGs 
(Farrah et al. 2001) is that ULIRGs are not a simple pre- 
cursor stage to optically selected QSOs, but instead are a 
diverse population whose evolution is driven by the progen- 
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Figure 2. AGN luminosity in the IR plotted against starburst luminosity in the IR. Triangles are AGN dominated systems, squares 
are Starburst dominated systems. Filled shapes are objects from the sample in this paper, outline shapes are taken from Table kL 



itor morphologies and the local environment. HLIRGs may 



de Zotti 199E). Gas and dust masses in HLIRGs have also 



simply be the high luminosity 'tail' of the ULIRG popula- 
t ion. Also worth no ting is that there is plausible evidence 
(Farrah et al. 2002) that HLIRGs may lie in very diverse 



been the subject of several previous studies, with derived 
dust masses lying in the range 10 8 < MduatjMp,) < 10 9 
arvainis, Antonucci fc Coleman 1992; |Powncs, Solomon~ 



environ ments, from isolated sources to Abcll class >2 clus- Radford 1995 ; Fraycr et al. 1998; |Frayer et al. 199£ ; Verma 



ters. Si ch a diversity of environment would give rise to a 
more diverse range of evolutionary paths. 

A recent model for the evolu tion o f HLIRGs in partic- 
ular is that of Rowan- Robinson (2000). According to this 
model, HLIRGs are ideal candidates for being 'primeval', 
or very young galaxies, as opposed to mergers between gas 
rich galaxies. Evidence in support of this is that HLIRGs 
all have very high star formation rates and a higher gas 
fraction than typically found in spiral galaxies. The data in 
this paper, coupled with previous results, allows us to ex- 
amine whether HLIRGs are very luminous galaxy mergers, 
or young active galaxies. 

The amounts of gas and dust in active galaxies have 
been examined by several authors and can be a useful di- 
agnostic. In the local Univer se, observations of nea rby ra- 
dio galaxies in the sub- mm (Knapp & Patten 1991) found 
some objects to have comparable dust masses (~ 10' M Q ) to 
those found in spiral galaxies, but also found other objects 
with dust masses as low as ~ 1O 4 M0. At higher redshifts, 
s tudies of samples of radio galaxies and radio quiet quasars 



( Hughe[j, Dunlop fc Rawlings 1997j; |McMahon et al~ 



1999 



Granato 1999[ ) derive dust masses 



Andreani, Franceschini 

in the range 10' < Mdust{MQ) < 10 a , with no appreciable 
dependence on redshift out to z ~ 5. These results imply an 
early accelerated phase of obscured star formation in high 
redshift primeval galaxies ( Franceschini et al. 1994 ; Mazzei 



et al. 2002) 



The dust masses derived for the sources in our sample 
are presented in Table ^. Assuming that the gas to dust 
ratio in HLIRGs is simi lar to that in ULIRGs, adoptin g a 
gas to dust ratio of 540 ( [Sanders, Scoville fc Soifer 1991 ) we 
derive gas masses in the range 10 lu -° < Mgas(Mo) < TO 115 . 
These gas and dust masses are comparable to those found 
by previous authors in other HLIRGs, and in other classes 
of active galaxy. A young, active galaxy may be expected to 
have both a very high rate of star formation (> 5OOM0yr _1 ) 
in order to form the bulk of its' stars in ~ lGyr, and a 
large gas reservoir {M gas ~ 1O 11 M0) to maintain this rate 
of star formation. Both these criteria are satisfied by all 
the objects in our sample. In addition, several of the host 
galaxies in our sample show both an elliptical morphology 
and no signatures of ongoing inte raction found in m ergers 
between equal mass spiral galaxies (Farrah et al. 2002). They 



are therefore excellent candidates for young, active galaxies. 
Some of our sample however have optical morphologies 



that do not fit with this scheme. HST imaging (Farrah et al 



2002 ) has shown that some HLIRGs, such as F00235+1024 
and F10026+4949 lie in interacting systems with similar 
morphologies to local ULIRGs. The most natural interpre- 
tation of these sources is that they are the high luminos- 
ity 'tail' of the ULIRG population. It is also conceivable 
that the starburst and AGN activity in the dynamically re- 
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Figure 3. Fraction of the IR luminosity due to AGN activity against the total IR luminosity. Filled squares are objects from this sample, 
outline squares are the additional HLIRGs from Table [l]. 



laxed systems in our sample are also triggered by galaxy 
interactions. The dynamical relaxation time of a merger be- 
t ween galaxies of comparabl e size is of the order ~ 10 9 years 
( Barnes fc Hernquist 1996 ; Dubinski, Mihos fc Hcrnquist 



1999 |) , idiereas the lifetime of a QSO (P.ees 1984; |Martini & 



Weinbe rg 200C ) and a starburst ( Genzcl ct al. 1998); [Thorn 



ley et a . 2000 ) is thought to be no more than 10 years. It is 



therefore just possible that a merger driven coeval starburst 
and AGN could be observed in a dynamically relaxed host, 
if the starburst and AGN activity were not triggered at or 
near the start of the merger. More likely though is that the 
objects in question are very young galaxies going through 
their maximal star forming period. We therefore conclude 
that the HLIRG population is comprised of both the most 
luminous end of the merger driven ULIRG population, and 
of a population of young active galaxies going through their 
maximal star formation periods. 



5.3 The Cosmic X-ray and Sub-mm Backgrounds 

The first three HLIRGs to be discovered were all shown to 



harbour obscured AGN (Hines et al. 1995). Just as Seyfert 2 
galaxies have been suggested to be Seyfert 1 galaxies viewed 
along a line of sight where the central regions are obscured 
by dusty tori, HLI RGs have bee n proposed as 'misdirected' 
QSOs, or QSO-2's (Hines 1998). A similar picture has also 



been prop osed for radio galaxies (Alighieri, Cimatti & Fos- 



bury 1994). Under this picture, HLIRGs would be classified 



as Hyp srluminous quasars if viewed along a polar line of 
sight, and as Hyperluminous Galaxies if viewed under any 
other orientation. 



We have shown that all the sources in our sample re- 
quire extremely luminous AGN to explain the total IR emis- 
sion. Two objects, IRAS 18216+6418 and ELAIS J1640+41 
were previously known to ha rbour AGN, since they are both 



et al. 



luminous in the hard X-ray (Pravdo & Marshall 1984; Morel 



2001). On the other hand, deep X-ray observation s of 
5 other HLIRGs (|Fabian ct al. 199$ jwilman et al. lgggh, in- 



cluding one (F00235+1024) from our sample, did not detect 
any of the objects, from which it was concluded that these 
objects were either atypically weak in the X-ray or that the 
central regions were obscured by columns Nh > 10 23 cm -2 . 
From the dust masses and AGN IR luminosities presented in 
this paper the latter conclusion is preferred. As the sample 
selection for our study is independent of AGN content, it 
seems likely therefore that many HLIRGs as a class contain 
a luminous AGN that is heavily obscured and can only be 
detected in the IR. 

The presence of highly obscured AGN in a large frac- 
tion of the most luminous /./MS galaxies has interesting im- 
plications for interpreting observations of the cosmic X-ray 
and IR backgrounds. Studies have shown that only ~ 10% of 
high-z bright sub- mm sources show signs of AGN activity a t 



other wavelengths (Barger et al. 2001; Almaini et al. 2001) 



Recent surveys have shown tha t the hard X-ray backg round 
is due almost entirely to AGN (Mushotzky et al. 2000), and 



there exist models for the X-ray background in which low 
luminosity AGN produce a hard X-ray spectrum via photo- 
electric absorption from an ob scuring starburst (Fabian et 



al. 1998 



Gunn fc Shanks 199S). Estimates of the AGN con- 



t ribution to deep sub-mm survey s and the IR background 
( Almaini, Lawrence fc Boyle 1998 ) assume that the fractions 
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Table 9. Predicted HLIRG SCUBA fluxes extrapolated to higher 
redshifts 



Name 


z= 


2.0 


z= 


3.0 




5850 


•S450 


•5850 


•S450 


F00235+1024 


5.3 


28.0 


6.7 


21.2 


07380-2342 


13.1 


59.4 


15.2 


47.5 


F12509+3122 


3.0 


15.3 


3.7 


14.9 


14026+4341 


2.9 


15.0 


3.5 


13.4 


F14218+3845 


5.6 


37.0 


7.2 


38.5 


18216+6418 


8.0 


40.5 


9.7 


34.7 



All fluxes are given in mjy. Predicted fluxes are given only for 
those 6 sources with constrained SED shapes. 

of obscured and unobscured AGN must fit with the observed 
cosmic X-ray background, and from this estimate that be- 
tween 10% and 20% of the 850/im SCUBA sources at lmjy, 
and a similar fraction of the IR/sub-mm background, are 
due to AGN. 

From our results however it is likely that, at least for 
the rest-frame most IR luminous sub-mm sources, this pic- 
ture is not true, and that there exists a population of heavily 
obscured AGN in sources with L ir > 10 12 ' 5 L Q . These AGN 
will contribute negligibly to the cosmic X-ray background. 
From the SEDs presented in this paper it is apparent that 
these AGN may supply a significant fraction, possibly up 
to 50%, of the contribution to the cosmic background ra- 
diation at \ bs > 300/^m from the most luminous sub-mm 
sources. Worth noting however is that these SEDs also pre- 
dict that the 850/xm background from the most luminous 
sub-mm sources will in all cases be due to star formation, 
with a negligible AGN contribution. 



5.4 Implications for Sub-mm Surveys 

5.4-1 Coeval starburst and AGN activity 

There exists substantial observational evidence that, at both 
low and high redshifts, the main episode(s) of star forma- 
tion in galaxies are linked with QSO activity. At low redshift 
there is an observed correlation between black hole mass 



and bulge mass ( Magorrian et al. 1995 ; Kormendy & Geb 



hart 2001) with a linear scaling between the two quantities 
( [McLure fc Dunlop 2002} ). At high redshift studies of emis- 
sion and absorption lines in high-z QSOs typically show that 
they reside in metal rich (at least IZq) environments, and 
t hat the more luminous objects may be more metal rich 
(Hamann & Ferland 1999). Sub-mm continuum observations 



of high-z QS Os generally show the presence of large m asses 



of cold dust ( pmont et al. 1996| ; |Benford et al. 1999[ ). This 
emission from cold dust is generally ascribed to st ar forma- 
tion, alt hough QSO activity is also a possibility (Carilli et 



al. 200C). The high metallicities in these systems is often 
seen to indicate that the optical QSO occurs after the bulk 
of star formation has occ urred in the system. A recent study 
of QSOs and their hosts (Granato et al. 2001 ) proposes that, 



in an optical QSO, the star formation is effectively at an end 
as the QSO power and feedback fro m supernovae com bine 
to blow out any nuclear starburst (silk & Rees 1998). An 
optical QSO activity and high star formation rates are thus 
predicted not to be coeval. A similar study ( Archibald et 



al. 2001) predicts that, unless a massive 'seed' black hole is 
already present, luminous AGN activity cannot arise until 
0.5Gyr after the first stars form. An independent estimate 
for the formation of a > 10 s Mq black hole in non-nucleated 



sources (Taniguchi, Ikeuchi & Shioya 1999) gives a timescale 
of lGyr. Therefore, this predicts that there is a time lag be- 
tween the peak in sub-mm luminosity (due to the starburst) 
and mid-IR luminosity (due to the AGN) of ~ 10 8 years. 
AGN activity is predicted to only occur in the host once 
~ 90% of the stars in the host have formed and the sub-mm 
luminosity is only ~ 2 5% o f its peak value. The predictions 
from Archibald et al. ( 2001) are that AGN in luminous high- 
z sub-mm sources are ~ 1000 times less luminous than the 
most luminous AGN found in other systems. 

The results from our study however suggest that, at 
least for the most luminous IR galaxies, coeval AGN activ- 
ity and massive starburst activity exist in many sources, 
including optical QSOs. Three objects in our sample 
(F12509+3122, F14218+3845 & ELAIS J1640+41) are opti- 
cal QSOs containing an extremely luminous AGN together 
with high rates (> 5OOM0yr _1 ) of star formation. The com- 
puted star formation rates are so high that it is unphysical 
to argue that they are not at or near the peak rates of star 
formation in these systems. In addition the strong corre- 
lation between the starburst and AGN luminosities in our 
sample, implying a common factor governing the luminosi- 
ties, suggests that the starburst and AGN phases cannot be 
widely separated in time in HLIRGs. 

5.4-2 The nature 0} high-z sub-mm sources 

Two SC UBA galaxies have been previously shown to be 
HLIRGs Jlvison et al. 2000| ). A sample of HLIRGs with se- 
lection unbiased toward AGN allows us to examine what 
fraction of HLIRGs as a class would be detected by SCUBA 
surveys if they lay at higher redshifts. Analysis can only be 
done for those six sources in the sample (Figure |l|) with con- 
strained SED shapes. The fluxes predicted to be observed 
by SCUBA for these sources if the sources lay at z = 2 
and z = 3 are given in Table |^. At these redshifts all the 
sources would have observed frame 850^im and 450/im fluxes 
comparable to the range of 850/im and 450/im fluxes for 
sources found in recent sub-mm surv eys (|Bargcr ct al. 1991 



Eales et al. 2000 



Scott et al. 2002 



Fox et al. 2002). Thus 



it seems plausible that at least some fraction of the high-z: 
sources found in sub-mm surveys are similar in nature to the 
galaxy populations found in the HLIRG population, namely 
a mixture of ULIRG-like galaxy-galaxy mergers and young 
active galaxies. It is thus plausible from our results that 
many of the sources discovered in sub-mm surveys harbour 
both a dusty starburst and a heavily obscured AGN which 
contribute comparably to the total rest-frame IR emission. 
Our study has also shown that, at rest-frame wavelengths 
> 50/im the emission from our sample is starburst domi- 
nated. It is therefore a safe assumption that high-z surveys 
sampling rest-frame wavelengths longward of 50/xm are for 
all sources sampling the starburst dominated region of the 
SED. 

It is, however, important to note that X-ray obser- 
vations of HLIRGs and X-ray surveys of high-z sub-mm 
sources imply that the two populations are not identical, and 
that there has been evolution between the two epochs. Of the 
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seven HLIRGs at z < 1 that have deep X-ray observations 
two (E J1640+41 & 18216+6418) are detected, correspond- 
ing to a detection fraction of 28%. Conversely, dee p Chandra 
observations of the SCUBA 8mJy survey regions ( Almaini et 
al. 200 J) detect only 2/36 sub -mm sources (assuming a de- 
tection threshold of 3.5cr for the sub-mm sources), or ~ 6%. 
When combined with our results this gives rise to two possi- 
bilities. The first is that the fraction of sources that contain 
luminous AGN is comparable at z ~ 1 and z ~ 3, with a 
greater mean obscuration level at z ~ 3. The second possi- 
bility is that the fraction of sources that contain luminous 
AGN has declined from z ~ 3 to z ~ 1. Further deep obser- 
vations of HLIRGs and high-z sub-mm sources in the X-ray 
and in the IR will be required to differentiate between these 
two possibilities. 



6 SUMMARY 

We have presented sub-mm photometry for 11 Hyperlumi- 
nous Infrared Galaxies selected purely on the basis of their 
IR emission, and used radiative transfer codes for starbursts 
and AGN in conjunction with IR photometry from the lit- 
erature to examine the power source behind the IR emission 
from 10 of the objects. Our conclusions are: 

(1) In all of the sources both a starburst and AGN com- 
ponent are required to explain the total IR emission. The 
mean starburst fraction is 35%, with a wide range span- 
ning 80% starburst dominated to 80% AGN dominated. In 
all cases the starburst dominates at rest-frame wavelengths 
longwards of 50^im, with associated star formation rates of 
at least 5OOM yr _1 . 

(2) The trend of increasing AGN fraction with increasing IR 
luminosity observed in the ULIRG population is observed to 
peak in the HLIRG population. The correlation between the 
starburst and AGN IR luminosities in the AGN dominated 
systems implies that the luminosities of both starbursts and 
AGN in HLIRGs are governed by common physical factors, 
the most plausible of which is the available quantities of gas 
and dust in the system 

(3) Comparisons between the fractional AGN and starburst 
luminosities, derived dust masses, and previously published 
HST imaging suggests that the HLIRG population is com- 
prised both of mergers between gas rich galaxies, as found in 
the ULIRG population, and of young active galaxies going 
through their maximal star formation periods whilst har- 
bouring an AGN. 

(4) The presence of a luminous AGN in all of the ob- 
jects in our sample, coupled with the non-detection in the 
X-ray of most HLIRGs, implies that many galaxies with 
Li r > 10 12 ' 5 Lq harbour a heavily obscured AGN. Although 
the contribution from these bright sub-mm galaxies to the 
total cosmic sub-mm background is small, the contribution 
from heavily obscured AGN in the brightest sub-mm sources 
to the cosmic background radiation at \ b 3 > 300/im may 
therefore be higher than estimates based on constraints from 
the cosmic X-ray background. The 850^m background from 
the most luminous sub-mm sources will however in all cases 
be dominated by star formation. 

(5) The detection of coeval AGN and starburst activity in 
our sources coupled with the strong correlation between 
starburst and AGN luminosity implies that, for the most lu- 



minous sub-mm sources, star formation does not cease once 
the source harbours a luminous AGN. We infer that star- 
burst and AGN activity, and the peak starburst and AGN 
luminosities, can be coeval in active galaxies generally. 
(6) When extrapolated to high-z our sources have compara- 
ble observed frame sub-mm fluxes to objects found in sub- 
mm surveys. At least some high-z sub-mm survey sources 
are therefore likely to be comprised of similar galaxy pop- 
ulations to those found in the HLIRG population, namely 
galaxy-galaxy mergers and young active galaxies. It is also 
plausible from these results that high-z sub-mm sources har- 
bour heavily obscured AGN. It should be noted that ~ 28% 
of HLIRGs at z < 1 are detected in the X-ray, whereas only 
~ 6% of high-z sub-mm sources are detected in the X-ray. 
This implies some evolution between the two epochs. Either 
the mean AGN obscuration level is greater at z ~ 3 than at 
z ~ 1, or the fraction of IR- luminous sources at z ~ 3 that 
contain AGN is smaller than that at z ~ 1. 
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